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Abstract 

Several shell deformation models are developed for use in crashworthiness analysis of rotationally 
symmetric structures. These models use analytical techniques to predict the crushing force 
versus axial crush distance characteristics of both a rigid- plastic, hemispherical shell and an 
elastic, cylindrical shell loaded axially by a rigid flat plate. Additional methods are proposed to 
determine the effects of cutout sections and internal stiffening members on the crushing force 
capacity of the shells. These proposed methods are applied to determine the energy absorption 
capability of the composite/metal nose structure of a mini submarine subjected to a head-on 
impact with a flat rigid wall. The nose structure is composed of a stiffened composite shell that 
covers and is attached to the metal, forward, hemispherical portion of the pressure hull. The 
complex structure is simplified to a rotationally symmetric shell model and modified to account 
for the effects of stiffening elements and cutout sections. Laminated plate theory, a progressive 
composite failure method, and the models developed are then used to determine the structure's 
energy absorption capacity. The fairing structure is determined to be capable of absorbing the 
kinetic energy associated with an impact from an initial vehicle speed of 7.5 knots. Further, the 
entire nose structure (the fairing and forward pressure hull) is predicted to be able to absorb 
the kinetic energy from a 26.4 knot vehicle impact. 

Thesis Supervisor: Tomasz Wierzbicki 
Title: Professor of Applied Mechanics 

Thesis Supervisor: James H. Williams. Jr. 

Title: SEPTE Professor of Engineering 
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Chapter 1 



Introduction 



The Advanced Seal Delivery System (ASDS) (Figure (1-1)) is a highly sophisticated, mini- 
submarine 19.8 m long that displaces 60 metric tons submerged. It is designed to transport a 
team of Navy SEALs from a host submarine to a tactical point of insertion and back. While 
the transported SEALs are a formidable offense weapon, the ASDS itself has no offensive or 
defensive capabilities and must rely on stealth and its ability to ’’hide' in the water column 
to ensure the safe transport of its personnel. As such, the ASDS will be operating in close 
proximity to the ocean bottom in order to mask its underwater signature. This near bottom 
operating strategy presents some problems, however. The most pressing of these is the increased 
risk of impact with submerged objects or structures. For the safety of the embarked crew, it is 
paramount that a determination be made as to the crashworthiness of the ASDS and whether 
operational restrictions need to be imposed. In effect, a safety envelope needs to be developed 
which takes into account the depth, bottom type and impact angle to give a maximum safe 
operating speed. A natural starting point for this analysis is to look at how a front end impact 
will affect the forward structure of the vessel. This analysis, however, presents no small task. 
The complexity of this problem quickly becomes apparent when the structure types and problem 
variables are viewed: 

• The structure of the nose section consists of various composite materials and steel in- . 



eluding: glass woven cloth, carbon cloth, kevlar. vinylester, rubber, syntactic foam, and 
HY-80 steel. 



Tail Section Aft Bod y 




Figure 1-1: External Isometric View of the ASDS 



• The variables of the problem include: depth, impact angle, speed, and bottom type. 

Additionally, the interaction of the structural material and modes of failure during defor- 
mation and fracture are not well understood or predictable to a high degree of accuracy. In 
light of these complexities, using model testing to develop some empirical means of relating 
the structural behavior during a collision to key parameters such as impact angle, velocity, or 
crush distance would be the optimal analysis tool. However, due to the cost associated with 
such tests, this method is not feasible. Moving to the other end of the spectrum, a computer 
simulated "crash” test of the ASDS could, in theory, be used to ascertain the desired safety 
envelope at a substantially cheaper price and in an accelerated time frame. While modeling 
the full ASDS via computer simulation seems plausible on the surface, a little digging reveals 
that even the most sophisticated computer tools have not evolved far enough to accurately 
predict impact responses for such complex, composite/metal structures. With this in mind, it 
is apparent that neither model testing nor computer simulation alone can solve this problem. 
A middle road between these methods may be the required path. In effect, a combination of 
structural simplification and computer and analytical modeling will be used to complete the 
analysis. With all of this said, the overall objective of the project will be to determine a safe 
impact speed for the vehicle such that during an impact the pressure hull will not rupture and 
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the crew will be able to survive the collision. 

As a first step in solving the problem, this thesis will analyze the energy absorbing potential 
of the composite nose fairing and forward pressure hull hemi-head structures during a head-on 
impact with a rigid wall structure. While this impact scenario is idealized and likely never 
to occur during actual operations, it will serve as a worst case impact and best case energy 
absorption scenario. This will be a worst case impact from the stand point that all kinetic energy 
will be absorbed by the structure in the form of deformation. No energy will be absorbed by 
the impacted structure or be maintained as vehicle kinetic energy. This is a best case energy 
absorption scenario in that the impact will allow the full utilization of the bow structure to 
absorb energy before the pressure hull is impacted. The reasoning behind looking at such a 
case first, is to ascertain the maximum energy absorption potential of the composite structure. 

This analysis will progress in three phases: 

• Deformation model development. 

• Structural modeling of the ASDS nose structure. 

• Energy absorption analysis. 

The models developed include deformation relationships for a rigid-plastic, hemispherical 
and an elastic, cylindrical shell loaded axially by a rigid flat plate. A method is also included to 
approximate the reduced load carrying capacity caused by cutouts in the shell structure. Next, 
the composite shell of the nose fairing will be modeled using simplified, rotationally symmetric 
shells (i.e.- cylinders, frustum, and hemispheres). Internal stiffening members are incorporated 
into the model by “smearing" the effective area of each stiffener across portions of the inner 
surface of the shell to increase the effective thickness of the structure. Cutouts in the shell for 
components, such as thruster doors, video equipment, and lights are also modeled. Finally, the 
deformation relationships developed will be applied to the corresponding sections of the model 
to produce a force verse displacement relationship of the structure as it is crushed axially by a 
rigid flat plate. The energy absorption capability of the entire structure will then be determined 
and translated into the vehicle speed needed to produce such an impact. 
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Chapter 2 



Development and Discussion of 
Analysis Tools 



2.1 Crushing Response of a Hemispherical Shell Between Rigid 
Plates 

2.1.1 Prior Research 

Much attention has been devoted to the study of the energy absorbing characteristics of thin- 
walled. metal, dome structures [l]-[7] . The problem of large deformations of rigid-plastic, spher- 
ical shells being compressed between rigid plates was first studied by Updike [1]. A simplified 
yield condition and geometry was used to obtain a sequence of limit loads, thus producing a 
force-displacement relationship valid from a displacement of a few shell thicknesses to about 
l/10th the shell radius. Wierzbicki [4] used similar geometric and deformation assumptions 
to derive a force-deflection relationship based on the optimization of bending and membrane 
energies valid theoretically until the plate deflection reaches half the radius of the shell. Com- 
parison of [1] and [4] show that the force-displacement results predicted by each differ only 
slightly. Experimental testing of hemispherical shells was reported by Kitching [5] and Kin- 
head [5]. Kitching et al studied the quasi-static loading of hemispherical shells and observed 
experimentally three different stages of deformation for shells with radius to thickness ratios 
ranging from 36 to 460. Later, Kinkead performed experimental testing on shells with radius 
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Figure 2-1: Geometry of Hemispherical Deformation Model. 

to thickness values ranging from 8 to 32 and identified two vice three deformation stages. The 
loading characteristic observed by Kinhead were compared to values predicted by the methods 
of [1] and [4] and large differences were found. Gupta [7] developed another prediction method 
using a rigid plastic analysis which relies on two calibration constants that must be determined 
from curve fits of experimental data. The method proposed by Gupta showed good correlation 
with experimental results. This analysis develops a means of predicting the load-compression 
relationship of hemispherical shells using a rigid plastic analysis technique with an emphasis 
on creating a simple, analytical model whose accuracy is improved over that proposed in [1] 
and [4]. The model developed will utilize the same deformation model used by [1] and [4], but 
calculate the membrane energy using a simple strain averaging method. Additionally, the re- 
sulting relationship will be capable of analyzing material with differing compressive and tensile 
stresses. 



2.1.2 Deformation Process 
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The spherical shell of radius R will be loaded from the top by a moving rigid plate with 
the base of the shell pinned to a fixed rigid plate. As the shell is loaded axially, a spherical 
dimple forms which is an inversion of the original shell curvature as shown in Figure (2-1). In 
this deformed state, surface AB remains rigid and undeformed. At point B, a circumferential 
plastic hinge increases the curvature of the shell, followed by rigid body translation between B 
and D and a reversal of curvature by the plastic hinge at D. The shell is subsequently unloaded 
and rigid body translation results on surface DE. As the shell deforms, the hinge circles at B 
and D move progressively outward. Deformation continues in this manner until the inverted 
portion of the shell at point E contacts the fixed rigid plate. 

2.1.3 Deformation Model 
Geometry 

This analysis begins with a look at the various geometrical relationships that can be attained 

from the deformation model in Figure (2-1). First, the radius to the center of the toroidal 

surface (a) and the displacement of the moving rigid plate (<5) can be expressed in terms of the 
shell radius ( R ) and deformation angle (o) as 

a = R sin 6, and (2.1a) 

6 = H(1 — cos<p) (2.1b) 

By differentiating Equations (2.1a) and (2.1b), the radial and axial velocities of the toroidal 
section and rigid plate become 



da 

dt 


do 

= R(cos6)—r- , and 
dt 


(2.2) 


d6 

dt 


o / • \ do 

= R(sind) — 


(2.3) 



By approximating cos (p(t) using the first two terms of its Taylor expansion in Equation (2.1b) 
and rearranging, a simple equation relating the deformation angle to the rigid plate displacement 
can be deduced. 

(2.4) 
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Turning now to the toroidal surface BCD , the average plastic hinge length (L) and toroidal 
surface area (.4) can be calculated by the Equations (2.5a) and (2.5b). 

L = 27ra, and (2.5a) 

A = 4tt rao (2.5b) 

Additionally, the radial and rotational velocities of the plastic hinge at point B can be taken 
as 



V 

d9 

dt 



do 

{ ~di 



V _ Rdo 
r r dt 



(2.6) 

(2.7) 



Finally, in order to simplify subsequent calculation, surface BCD will be approximated as a 
parabolic surface of the form 



y 



1 £ x 2 
2 r sin ip 



+ r(l — cosd) 



with the x and y axes oriented as shown in Figure (2.2.3). 



(2.8) 



Derivation of Membrane Stresses and Fully Plastic Bending Moment 

Next, the membrane force (N) and fully plastic bending moment (A/) relationships for 
a material having differing compressive ( a c ) and tensile (at) flow stresses will be developed. 
As shown in Figure (2-2), when a c ± a t the neutral axis will shift away from the geometrical 
center of the cross section by an amount 77 . By using the fact that the neutral axis of a specimen 
subjected to a bending moment is located where the membrane force (AT) is equal to zero, we 
are able to calculate a general expression for 77 in terms of shell thickness ( t ). a c . and a t as 
shown below 

fi-v , , 

N = / a t dz — / a c dz = 0 (2.9) 

Jo Jo 

where both at and a c are taken to be positive.. Upon integrating and solving for 77 , Equation 
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Figure 2-2: Membrane Stress and Fully Plastic Bending Moment Representation for Material 
having Differing Compressive and Tensile Flow Stresses. 



(2.9) becomes 



1 <r c -<Tt 
V = o*: 



2 cr c + (7 t 

With 77 calculated, general expressions for N and M can be derived as follows: 



( 2 . 10 ) 



.V = f 
M = J 



(j c dz 

h+n 



a t zdz ■ 



-f 



(JrZdz 



Note that the regions of membrane deformation are compressive in nature while the bending 
forces are both tensile and compressive. After integrating, substituting Equation (2.10) for 77 . 
and setting a = 2 a °+ at ■ the total membrane force and fully plastic bending moments become 



N 



M 



<7 c t 




&c + &t 




( 2 . 11 ) 

( 2 . 12 ) 



In order verify the validity of Equations (2.11) and (2.12). a c can be set equal to a t which 
results in a symmetric stress profile throughout the specimen. In such a case Equations (2.11) 
and ( 2 . 12 ) reduce to the standard expressions for the fully plastic membrane force and bending 
moment. 
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Energy Balance 



The principle of virtual velocities wall be used to relate the rate at which energy is being 
introduced to the system to the rate at which it is being absorbed. Equation (2.13) states that 
the rate at which energy is applied to the shell must exactly match the rate at which it is being 
absorbed by the bending and membrane effects in the shell. 



(1) (2) (3) 

p dS_dEb dE m 
dt dt dt 



(2.13) 



The remainder of this section will further develop each of the terms in Equation (2.13). Starting 
first with the left hand side of Equation (2.13). term (1) can be rewritten in a different form 
by using Equation (2.3) as follows 

pd i = PR ^ d i ^ 



The second term of Equation (2.13) represents the rate at which energy is absorbed by the shell 
due to bending processes. Using the deformation model shown in Figure (2-1), it is plainly seen 
that bending energy will be absorbed at two discrete locations - plastic hinges B and D. This 
component can be quantified by the following expression 



dE b 

dt 



E 

i=B,D 



(2.15a) 



Equation (2.15a) can be further simplified by making the following assertions 



• the magnitudes of the bending moment and rate of angular rotation at hinges B and D 
are equal. 



• the total hinge length 
length. 



Li can be rewritten as 2 L (where L is the average hinge 

i=B,D 
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Figure 2-3: Deformation of Surface BCD. 



Using these simplifications and Equations (2.1a). (2.5a). and (2.7), Equation (2.15a) can be 
expressed in the following form: 



dE b 

dt 



t 2 . ,Mdo 

4 ttR sin o — 

r dt 



(2.16) 



The final term in Equation (2.13) accounts for the membrane energy dissipated on surface BCD 
and can be written as: 



^ =NA £ 

dt dt 



(2.17) 



where ^ is the instantaneous rate of strain experienced by surface BCD. As a simplification, 
the average strain rate, will be used in place of the instantaneous strain rate. Figure 
(2-3) shows a detailed view of the deforming surface BCD. The average strain of this surface 
will be calculated by dividing the surface into two subsurfaces, calculating the average strain 
on each, and then averaging the strains from each region to get the final average strain. This 



statement can be expressed mathematically as 



Sav 



g/qv + £//av 
2 



(2.18) 



where £j av and c// au are the average strains experienced as surfaces BC 0 (region 1) and C 0 D 0 
(region 2) deform to BCand CD respectively. Looking first at region I. the strain can be written 
in terms of initial and final curvatures of the BC. 



si = e If - e Io 



1 

2 





In general, the average strain in region I can be formulated as: 



a, 



i r s 

sin 0 Jo 



Eidx 



(2.19a) 



( 2 . 20 ) 



Upon substitution of Equation (2.19a). integration and simplification. Equation (2.20) reduces 
to 

e /„, = 1<* 2 (2.21) 

The strain in region II requires an additional step to calculate. While the strain in region I 
is purely compressive for the entire deformation, the strain in region II is initially compressive 
for the first part of deformation and then tensile for the remainder of the deformation. So, 
the average strain will be calculated for the compressive phase as surface C 0 D 0 is deformed to 
an intermediate flat surface C*A and also for the tensile phase as the shell further deforms to 
surface CD. With this, the average strain in region II can be calculated in the same manner 
as in region I. The details of the calculation follow. 



c //at. = s i / a , 



Tipressiue 

rr sin0 






Sll av 



T Sin O jo [ compressive ^^°comprejjive ) I B tensile ^°tensWe ^ ] 

■ ^r[ G'-M-Ki) 



= 3 6 ' 



dx 



( 2 . 22 ) 



18 



With the average strain calculated in both regions, we can now determine the average strain 
of the entire surface BCD by substituting Equations (2.21) and (2.22) into Equation (2.18). 
Subsequently, the average strain rate can be calculated by simple differentiation of the average 
strain as follows. 



Eav 

dSav 

~df 



,do 

dt 



(2.23) 

(2.24) 



Now a final expression for the rate at which energy is absorbed by membrane forces can be 
realized by substituting Equations (2.5b). (2.1a) and (2.24) into Equation (2.17). 

dE™ — 47rr NR sin oo 2 (2.25) 
dt dt 



New expressions have now been formulated for each term in Equation (2.13). By substituting 
the new expressions (Equations (2.14), (2.16), and (2.25)) into Equation (2.13) and solving for 
the instantaneous crushing force (P), the following result is obtained 

(2,26) 

Further, the mean crushing force (P) can be calculated and minimized with respect to the 
toroidal radius (r) to obtain an optimal toroidal radius. 



P = 
dP 

dr 





(2.27) 



(2.28) 



By substituting the optimal toroidal radius back into Equation (2.27) and using the fact that 
at the optimal radius the membrane and bending contributions are equal, the mean crushing 
force can be expressed more simply as 
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Popt — 4 oR 




(2.29) 



Equation (2.29) can be rewritten in terms of plate displacement and shell thickness by using 
Equations (2.4), (2.11), and (2.12). 



2.2 Derivation of Cylindrical Tube Equation 

2.2.1 Prior Research 

The deformation and energy absorption characteristics of thin- walled, composite, cylindrical 
shells have been studies for many years by many researchers. [8]-[15j. Farley [10]-[13] and Hull 
[14] have performed extensive testing on composite shell samples composed of various mate- 
rials using different lay-up conditions. Both quasi-static and dynamic loading situations were 
investigated. These studies did much to further the knowledge base in the areas of determining 
crushing modes and their controlling mechanisms. They also showed that there was no explicit 
set of rules that could be generalized for all composites types. Material lay-up pattern, crushing 
speed, fiber volume fraction, and radius to laminate thickness ratio all play key roles in deter- 
mining the deformation characteristics of the shell. From these studies, complicated numerical 
models have been established for use in predicting the failure of very simple composite struc- 
tures. The degree of success realized by these models has varied in even simple situations. The 
vast majority of the work described has been directed toward understanding the deformation 
of composites on the micro- and macroscopic scales. Hanefi and Wierzbicki [15] took a more 
global approach to a similar deformation problem. They developed a simple model that was 
used to predict the energy absorption characteristics of a combined metallic/composite tube 
structure in axial compression. The model coupled a rigid plastic analysis of the metal tube 
with a compressive perfectly plastic and tensile elastic analysis of the composite material. The 




(2.30) 



Finally, the optimal instantaneous crushing force becomes 



Popt — 47 T\/ 2 cr c a 2 £ 2<52 



(2.31) 
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Figure 2-4: Deformation Model for Cylindrical Thin-walled Structure. 

predicted results matched well with experimental data. In this analysis, a global failure ap- 
proach will also be taken. A deformation model will be developed using elastic principles to 
determine the force deflection and energy deflection characteristics of a thin-walled, cylindrical 
shell loaded axially. Bending and in-plane stiffnesses will be step-wise varied to account for 
stiffness variations as plys fail within the structure. Finally, a simple result will be established 
for a material exhibiting constant and identical flexural and inplane stiffness properties. 

2.2.2 Deformation Process 

Figure 2-4 shows a representation of two buckling lobes of the cylinder during deformation 
where 6 denotes the axial crush distance, a is the radius of the cylinder, and 2 H is the half 
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length of the buckling lobe generated. As the cylinder is crushed axially by a distance <5, the 
cylinder will form an axisymmetric, elastic buckling pattern with a wave length equal to 2 H . 
Deformation takes place both in the form of bending and circumferential strain. The material 
will behave elastically until its elastic limits are reached with subsequent deformation occurring 
with the material in a progressively degraded state. The material degradation takes place 
through a step-wise reduction of the material's flexural and circumferential elastic moduli and 
material thickness. The process continues progressively through all buckling lobes formed. 

2.2.3 Deformation Model 

The principle of virtual work (Equation (2.32)) will be used as the starting relationship from 
which to derive an expression relating the crushing force ( P ) to the axial crush distance (<5). 
Equation (2.32) states that the energy applied to the cylinder through an external force (1) 
must be exactly equal to the energy absorbed by the structure in the form of circumferential 
(2) and bending (3) strains. 



(1) 


f ( 2 ) ( 3 ) 




PdS = 


J (crgdeg + a x de x )dV 


(2.32) 


(7Q = 


Eeg and a x = Ee x 


(2.33) 


PdS = 


J (Ecgdcg -r Ee x de x )dV 


(2.34) 



Because the deformation will be elastic, Equation (2.32) can be simplified using the elastic 
relationships between stress and strain (Equation (2.33)) to obtain Equation (2.34). Figure (2- 
5) shows a detailed view of the deformation of a single buckling lobe with parameters pertinent 
to the analysis labeled. Using this figure, expressions can be derived in terms of the variables 
shown for the incremental volume (dV), strains ( e x and eg). and incremental strains ( de x and 
deg). First, inextensibility will be assumed for the shell as it buckles, making dV constant over 
the deformation process and having the following value: 

dV = 4 ? ratH (2.35) 



99 



3 



Figure 2-5: Detailed View of the Deformed Cylinder Section. 
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Next, the circumferential strain can be expressed as the ratio of the outward deflection of the 
shell (w) to the original shell radius. The actual deflection is measured from the original surface 
of the shell to the deformed surface and varies from 0 at the being of the lobe to a maximum 
value of Wo at the midpoint, and then back to 0 at the end of the lobe. In order to simplify 
future calculations, a linear deflection profile {w(x) = jjii'o) is used to approximate the actual 
deflection profile as shown in Figure (2-5). Using this linear relationship, eg and deg can be 
expressed as 



eg 

deg 



w x 

— = ~T7 w o 
a aH 



aH 



dli'o 



(2.36) 

(2.37) 



The bending strain component can be determined by using the curvature of the shell section 
(k) and the position within the thickness of the material (z) as is shown in Equation (2.38). 
This variation of the bending strain throughout the shell thickness will be assumed linear with 
maximums occurring at | and and no strain occurring at the mid-plane. Further, plane 
sections will remain plane as shell curvature increases. 



z 




(2.38) 



However, the radius of curvature (R) is related to the half lobe length and the maximum 
circumferential deflection by Equation (2.39) 



R 2 -(R- Wo ) 2 

R 



H 2 -wl 
f r- 



2u.'o 



(2.39) 



Substitution of Equation (2.39) into Equation (2.38) yields the final expressions for e x and de, 
(Equations (2.40) and (2.41)). 



Now. Equations (2.35)-(2.37) and (2.40) can be substituted back into Equation (2.34) to obtain 
the following: 



PdS : 



( (Altai) 



a 2 H 2 



T 1 J E m (x) 2 dx + ^[(AitaH) J z 2 dzj w 0 dw 0 (2.42) 



If the flexural and membrane Young's moduli can be assumed equal and constant throughout 
the deformation process. Equation (2.42) can be further simplified to Equation (2.43 ). 




't_H 

a 



at 3 

H 3 



W 0 dlL' 0 



(2.43) 



Previous work [18] has proven that the elastic, buckling lobe length of a cylinder can be related 
to the shell radius and thickness by 



2H = 1.72 Vat 

H a Vat (2.44) 



Using Equation (2.44), Equation (??) can be reduced to 




(2.45) 



A final simplification can be made to Equation (2.45) by relating 6 to w Q . Using Figure (2-5) 
the following geometric relationships can be generated. 



sin(a) = ^ (2.46) 

cos(a) = 1 - (2.47) 
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By squaring and adding Equations (2.46) and (2.47). the relationships are combined in the 
following form 



After differentiating, rearranging and substituting Equation (2.44), the desired relationship is 
achieved 



This expression can now be substituted back into Equation (2.45) to obtain the final expression 
relating the crushing force applied to the axial deflection of the cylinder. 



As previously states, Equation (2.49) applies to materials with equal and constant flexural and 
circumferential Young's moduli. If these conditions are not met, the more general expression 
(Equation (2.42)) must be used. 

2.2.4 Correction for Open Sections 

Open sections or cutouts in a shell will effectively reduce the load required to deform the shell 
in the local region of the cutout. These cutouts could even change the local deformation mode 
of the structure or change the order in which the overall structure deforms. For simplicity, this 
analysis will include only the reduction in loading capacity for several reasons. The relationships 
previously developed. Equations (2.26) and (2.49), are valid only for the failure modes assumed: 
so, a change in failure modes (i.e.- from axisymmetric shell yielding to lobar buckling for ex- 
ample) will not be captured here in. While this seems like a fairly unreasonable assumption, 
the reader is reminded that the objective of this analysis is to determine the maximum energy 
absorption capability of the structure. The development of Equations (2.31) and (2.49) assume 
failure modes which produce this maximum. So. a shift to a different mode would produce a 
lower energy absorption. With the assumption that the failure mode remains unchanged at the 
cutouts, it does not make sense to consider the order in which parts of the structure collapse 
because the same amount of energy will be absorbed regardless of the collapse sequence. So. the 




(2.48) 





(2.49) 
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only effect that will be considered is the reduction of strength caused by the loss of material. 
This effect will be approximated by multiplying the applied load of an intact shell by a reduc- 
tion factor in the local vicinity of the cutout. This reduction factor will be proportional to the 
remaining shell area divided by the material present in an intact shell of the same dimensions 
(Equation (2.50)) 



F (6) corrected — F(.6) uncorrec t e d^{S) where 

r(A enioved material 



e(6) 



A,, 



•ompletc shell 



(«) 



(2.50) 



While this area ratio is part of the reduction effect, it does not capture all aspects. Other 
factors, such as cutout geometry, affect the load carrying capacity of the structure around 
the cutout. While these other effects are recognizable, a theoretical derivation of their effect is 
impractical. This effect must be determined experimentally. The experimental correction factor 
would be used to match experimental test data with theoretical predictions and is expressed 
as f in Equation (2.50). Since no experimental calibration was done during the course of this 
thesis. £ will be set equal to 1 until future evaluation determines its correct value. 



Hemispherical Dome 

For the hemispherical dome, the cutout sections will be assumed to be circular with diameter 
D as shown in Figure (2-6). From Figure (2-6) it is apparent that the correction factor can be 
stated in terms of the average plastic hinge radius (a(<5)) and angle o(6) as follows: 



(Average Hinge Diameter-Length of Open Section) \ 

Average Hinge Diameter / hemisphere 

/ 2»ra(g)-^)a(g) y, /. _ 

V 2ira(6) ) ? ' ,em “ p ' ,ere \ 2v 



£ hemisphere 



(2.51) 




The task now becomes to express o(6) in terms of known quantities. Figures (2-7) and (2-8) 
show more detailed views of the side and top of the shell respectively with pertinent dimensions 
labeled in each figure. 
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Figure 2-6: Top View of Hemispherical Dome with a Single Circular Cutout Section. 
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Figure 2-7: Geometrical Relationship of Parameters Used in the Derivation of the Correction 
Factor for Open Sections in a Dome Structure (side view). 

From these figures, the following geometric relationships can be established: 



Al = r Q + DcosO - a(S) 
Al 

A2 = 

cost/ 

A3 = D-A2 



D D 

A4 = — - A3 = A2 - - 

2 2 




(2.52) 



(p(6) can now be expressed in terms of the initial radius at which the cutout is encountered (r D ), 
cutout diameter (D), cutout inclination angle (#), and average plastic hinge radius as shown in 
Equation (2.53). 
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Figure 2-8: Geometrical Relationship of Parameters Used in the Derivation of the Correction 
Factor for Open Sections in a Dome Structure (top view). 
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o(<5) 

o(S) 



2 arcsin 



•m 

( / f r a +D cos 0—a(6)\ ( r\ { r a +D cos 0—a(6)\\ 

yl cose Jr~l )) 

a(S) 



(2.53) 



Next. a(<5) can be related directly to the radius of the hemisphere (R) and the crush distance 
of the shell (<5) by the following: 



a(6) = y/R 2 - (R - S ) 2 = yJb(2R - 6) (2.54) 

Finally, substitution of Equation (2.54) into (2.53) and the result into (2.51) yields the desired 
result for the correction factor (Equation (2.55)) 



e(S) = 



1 — — arcsin 



( r 0 +Dcase-y/6(2R-6)\ ( ^ ( r 0 +DcosB-y/6(2R-6) 

1 cos 6 /V V COs6 



s/H2R - 6) 



)\ 



£ hemisphere 

(2.55) 



Cylindrical Shell 

The same principles described in Section (2.2.4) will also be used to estimate the effect of 
cutout sections in the cylindrical section with the exception that the cutout shape will be much 
simpler. The cutouts in the cylindrical section will be rectangular in shape, and thus, the 
material removed will not vary with axial crush distance. So, using Equation (2.50) and again 
assuming that £ = 1. the correction factor of the cylinder can be expressed as 

2i:a — mpa no , . „ _ 

£ = 2^ = 1--^: where (2.56) 

n = Number of cutouts 

tp = Angular extent of cutout 
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A stiffener 





Figure 2-9: Stiffener Smearing Process. 



2.2.5 Effect of Shell Stiffeners 

The energy dissipation effects of shell stiffeners can be taken into account using Equations (2.31) 
and (2.49) by assuming that these structural members are "smeared'" across the surface of the 
shell so as to produce a new effective shell thickness. This process is illustrated in Figure (2-9). 
If the material properties of the stiffeners are identical to those of the shell material, their cross 
sectional area is distributed across the surface of the shell such that the new effective thickness 
( t e ff ) can be determined by the following expression 



the shell radius, and t s h e ii is the shell thickness. If the shell and stiffener materials have different 
properties, the force distribution over the effective cross section must be made consistent with 
the force distribution prior to the "smearing" process by using the following relationship 



n 




(2.57) 



where A t is the cross sectional area of the i th stiffener, n is the total number of stiffeners, a is 



f*eff — Pshell P stiffeners 



(2.58) 
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By substituting the appropriate relationships and assuming that P e fj will be calculated usinj 
the shell properties. Equation (2.58) can be written as 



2tt at e f fE s h e u£ — 2'nat s h e uE s heii£ ^ AiE t c 

1=1 

which can be solved to determine the effective shell thickness 

1 n 

te fl = t shel i + 2^L^ 

Z = 1 

When E\ = E 2 = ... = E n in Equation (2.59). Equation (2.57) is recovered. 



(2.59) 
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Chapter 3 



Modeling and Analysis of ASDS 
Composite Shell Structure 



With the analysis tools developed, focus now shifts to the modeling and analysis phase of 
the project. First, a model will be developed for the actual ASDS nose fairing structure and 
forward hemi-head of the pressure hull using simple, rotationallv symmetric, geometric shapes 
for the shell. Stiffeners and cutout sections will be incorporated using the techniques discussed 
in Section (2). Next, the analysis tools developed in Section (2) will be applied to determine 
the crushing force to axial crush distance relationship and energy absorption capability of the 
structure. Finally, a maximum impact speed will be calculated by equating the energy absorbed 
in the collision to the kinetic energy possessed by the vehicle just prior to the impact. 

3.1 ASDS Modeling 

3.1.1 Actual Structure 

The forward portion of the ASDS consists of a hemispherical pressure hull and a composite 
fairing structure. The hemispherical pressure hull section is constructed of HY-80 steel and 
is welded to the cylindrical portion of the hull to form a water tight boundary. The ASDS 
nose fairing (see Figure (3-1)) is a composite structure attached to the forward end of the 
pressure hull by eight titanium clevises. The structure has an overall length of 2.94 m and a 
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Figure 3-1: The Actual ASDS Nose Fairing Structure Showing Internal Components (Viewed 
from Aft Looking Forward). 
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Property 


A 


B 


C 


D 


E 


Ei (x 10 6 psi) 


3 


11.6 


4.7 


.0013 


16.4 


E2 (x 10 6 psi) 


2.7 


10.6 


4.7 


.0013 


16.4 


E3 (x 10 b psi) 


.55 


.37 


.45 


- 


- 


^1.2 


.15 


.124 


.15 


.499 


.31 


Fiber Content 


.5 


.5 


.5 


- 


- 


Tensile Strength 












Axial (X) (ksi) 


39.1 


75 


56 


- 


125 


Transverse(Y) (ksi) 


39.1 


75 


56 


- 


125 


Compressive Strength 












Axial (X # ) (ksi) 


33.8 


68 


14 


- 


125 


Transverse ( Y ) (ksi) 


33.8 


68 


14 


- 


125 


Shear Strength (S) (ksi) 


9.41 


7.5 


4.5 


- 


77 


A- E-Glass (1583)/Epoxy (7780) woven cloth 
B-Carbon woven cloth (6K 5HS satin weave IM7) 
C-Kevlar (K285)/ epoxy (7714) woven cloth 
D- Rubber core material 
E-Titanium devise material (Ti 6A1-4V) 



Table 3.1: Properties of Materials Used in ASDS Nose Fairng. 



pseudo-rectangular cross section with maximum dimensions of 2.21 m high and 2.06 m wide. 
The structure itself consists of an outer shell constructed from E-Glass (1583) /Epoxy (7780) 
preimpregnated (prepreg) woven cloth in a single mold process. Internally, the shell is stiffened 
by a network of longitudinal and transverse woven carbon cloth (6K 5HS Satin Weave IM7) 
beams which are bolted to the shell and interconnected by titanium clevis fittings. Attached to 
the front end of the shell is a sonar dome comprised of a BFGoodrich Rho-Cor material system. 
This system consists of Kevlar (K285)/Epoxy (7714) prepreg woven cloth with a rubber core 
material. The properties for the various material systems described are contained in Table 
(3.1). Cutouts have been made in shell and the sonar dome to accommodate such equipment 
as thrusters, the forward anchor, lights and video cameras. 



3.1.2 Simplified Model 

The actual structural design of the ASDS had to be simplified to rotationally symmetric 
shapes before the analysis tools developed in Section (2) could be applied. A representation 
of the model developed is shown in Figure (3-2). For the forward pressure hull, no modeling 
idealizations needed to be made since the structure was already rotationally symmetric and 
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9! 



Camera Cutout 



1.36 



5 .63 



28 DIA 



.91 R 
1.07 R 



.71 



Thruster Door Cutout 



1.13 



.99 R 



.98 R 



Anchor - 

Transverse Beams 

Forward Pressure Hull Hemi-head 



- Light Cutou 



Figure 3-2: Model of ASDS Forward Section (Dimensions in meters). 

hemispherical in shape. The pressure hull has a .99 m radius and an undisclosed thickness. The 
actual thickness of the structure is held as confidential by the U.S. Navy so as not to reveal 
the operational diving depth of the vessel. For this analysis, a notional thickness of 19 mm 
will be used. The sonar dome was modeled as a spherical dome section using the axial depth 
of the actual sonar dome and a base radius equal to the average of the height and width of 
the actual sonar dome. The thickness and lay-up construction was taken to be that of the 
actual structure. Three cutout sections simulating the locations of the camera and lights were 
located on the model 120° apart at the same radial distance. The actual and model sonar 
dome dimensions are summarized in Table (3.2). The shell structure was modeled using a 
combination of a frustum for the forward portion of the shell that transitions into the sonar 
dome and a stiffened cylinder for the after section. The radius of the cylindrical section was 
taken to be the average of the maximum height and width of the actual structure. The frustum 
was used as a transition section between the sonar dome and the cylindrical section and uses 
the corresponding radius at each end to fix its dimensions. For the frustum section, the shell 
thickness was taken to be that of the actual material. The thickness of the cylindrical section 
was taken to be an effective thickness calculated using Equation (2.59) to account for the added 
structural stiffness of the longitudinal carbon beams. The transverse beams were not utilized 
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Actual 


Parameter 


Model 


.64 


Axial Depth (m) 


.635 


1.91 


Height (m) 


- 


1.75 


Width (m) 


- 


1.83 


Average (m) 


- 


- 


Radius (m) 


.91 


9.6 


Thickness (mm) 


9.6 


Camera Cutout 


.28 


Diameter (m) 


.28 


0° 


Orientation 


0° 


.71 


Radial Distance (m) 


.71 


Light Cutout 1 


.28 


Diameter (m) 


.28 


-147° 


Orientation 


-120° 


.72 


Radial Distance 


.71 


Light Cutout 2 


.28 


Diameter (m) 


.28 


00 


Orientation 


120° 


.69 


Radial Distance 


.71 



Table 3.2: Comparison of Actual and Modeled Sonar Dome Dimensions. 

in the effective thickness calculation because their energy absorption in an orthogonal impact 
is would be minimal. Instead, these beams were used to define nodal points for buckling waves. 
Table (3.3) summarizes the actual and modeled dimensions of the shell section. Two cutout 
sections 1.13 m high by .81 m wide were included in the cylindrical section to account for the 
openings created by the forward thruster doors. (It is assumed that during a collision imminent 
situation, the thruster doors would be open and the thrusters in use.) Additionally, the titanium 
clevises (used to attach the internal beams together and to the pressure hull) and the anchor 
were treated as rigid objects. 



3.2 Analysis 

3.2.1 Impact sequence 

As stated previously in Section (1). this analysis will investigate an head-on impact with a. 
flat, rigid wall. Structural deformation will initiate at the bow as the sonar dome is impacted 
first. The deformation will progress .635 m through the sonar dome, .95 m through the frustum 
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Actual | Parameter | Model 


Cylindrical Section 


1.36 


Length (m) 


1.36 


2.21 


Height (m) 




2.06 


Width (m) 




2.14 


Average 




- 


Radius (m) 


1.07 


12.8 


Thickness (mm) 


35 


Frustral Section 


.95 


Length (m) 


.95 


- 


Fwd Radius (m) 
(Sonar Dome Radius) 


.91 


- 


Aft Radius (m) 
(Cylinder Section Radius) 


1.07 


10.1 


Thickness (mm) 


10.1 



Table 3.3: Comparison of Actual and Modeled Nose Shell Section Parameters. 



section of the nose shell, and .18 m into the cylindrical section before the anchor makes contact 
with the wall. At this point, some of the crushing load is transferred to the pressure hull via 
the anchor and simultaneous deformation of the cylindrical section and the pressure hull results 
for another .66 m. The deformation process ends when the after most transverse beams and 
titanium cleaves make contact with the wall. 



3.2.2 Sonar Dome 

Equation (2.31) will be used to determine the crushing force required to deform the structure 
throughout the progressive crushing of the sonar dome. Implicit in the use of this equation is 
that the structure can be modeled using rigid-plastic techniques. While it is understood that in 
general composite structures behave in a more brittle than ductile way, kevlar composites have 
been shown to buckle under axial loading in modes similar to those seen in metal structures 
[19] . This behavior coupled with the fact that the structure has a rubber core was used to 
justify the use of this analysis technique on this component. 

The sonar dome is composed of 22 plys of kevlar cloth laid-up symmetrically sandwiching a 
3.81 mm rubber core ([90. 0]s, 0, rubber core. 0, [90, 0 ) 5 ) for a total laminate thickness of 9.4 mm. 
The compressive and tensile flow stresses were calculated using a combination of laminated 
plate theory and a progressive failure analysis [21]. The detailed calculations are shown in 
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Appendix (A) and, for brevity sake, only an summary of the process and the result will be 
included here. First, the material properties of the ply (listed in Table (3.1)) were used to 
determine the undamaged 0° ply stiffness characteristics. These properties were then modified 
using techniques of [21] to determine the reduced stiffness of the 0° ply with matrix and fiber 
damage. The stiffnesses were calculated using the following relationship: 






— &L-. 
1— I'it'2 

E o 

\—V\ U 2 



^E 1 



...E 2 
1 — 1/2 



0 

0 



0 0 e 3 



(3.1) 



The stiffnesses were then rotated to a 90° ply orientation. The stiffness of the rubber core 
material was also determined from Equation (3.1) by setting E\ = Eo. rq = v 2. and £3 = 

The overall laminate stiffness was then calculated using the following: 



A* 

Ax 

h 

Q 



1 

^laminate 



(Ao°ho°Qo° ^900^900^90° “f“ A rubber hrubberQ rubber} where 



Total number of "x' ? plys 
Thickness of individual ply 
Stiffness of individual ply 



Next, stress space strength parameters were calculated from the material test data provided in 
Table (3.1) using Equations (3.2) and (3.3) and the generalized von Mises value of — ^ for F± 2 - 
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(3.2) 



(3.3) 



The strain space strength parameters and failure envelopes could then be determined through 
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the relationships in Equations (3.4)-(3.6). 



G u 


— P i,jQi,kQj.i 


(3.4) 


Gj 


= FrQij 


(3.5) 


GtjSiSj + GjS t 


= 1 


(3.6) 



A progressive failure model [22] was then used to approximate the stress strain relationship of 
the material under pure compressive and tensile loadings (oriented in the 0° ply direction). The 
progressive failure model allows for a rough prediction of post first ply failure behavior of the 
laminate. The model works by calculating the strain of first ply failure for the assumed load 
condition and then determines the stress using the undamaged laminate stiffness ({cr} = A*{c}). 
After first ply failure, the limiting ply is failed to a degraded state, a new laminate stiffness is 
calculated, and the next ply failure strain and stress is determined. A ply can be damaged at 
most twice - once by degraded matrix and the final by fiber failure. In certain cases, based on 
the properties of the plvs used, the degraded ply failure strains can be less than that of the 
undamaged ply. For this condition, the ply failure occurs after its undamaged failure envelope 
is surpassed. The tensile and compressive stress versus, strain estimates are shown in Figure 
(3-3). 

The compressive ( <7 ^ ) and tensile (c 7 to ) flow stresses were then determined to be 8.41 MPa 
and 53.35 MPa respectively using Equation (3.7). 

a 0 = — / crds (3.7) 

£ f Jo 

Now. the crushing force (P so »ar dome) and cumulative energy absorption (E SOUiir dome) of the 
sonar dome can be determined from Equation (2.31) and (3.8) with the cutout correction factor 
(Equation (2.55)) applied to the crushing force over the appropriate regions of the shell. 

E = f* P(6)dS (3.8) 

Jo 

The resulting crushing force and cumulative energy absorption profiles for an axial crush dis- 
tance (6) are shown in Figure (3-4). Here it is seen that the total energy absorbed for a crush 
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Figure 3-3: Compressive and Tensile Stress vs Strain Estimates Calculated using the Progressive 
Failure Method. 

distance of .635 m is .058 MJ. This figure also clearly shows the effect of the cutouts on the 
crushing force curve. 

3.2.3 Frustum Section of Nose Shell 

The frustum section is composed of a 22 ply. E- Glass (1583)/Epoxv (7780), woven cloth com- 
posite system laid up symmetrically ([90,45, 0, —45)2,90,45, 0,0, 45,90, [90, 45, 0, — 45]2) for a 
laminate thickness of 10.1 mm. The analysis of the frustum section was done by approximating 
the frustum as a cylinder with a radius equal to the average of the forward and after radii of the 
frustum and using the analysis discussed in Section (2.2.3). The details of these calculations are 
contained in Appendix (B). In order to use this technique, the appropriate circumferential (Eg) 
and flexural (E/) Young's moduli needed to be determined. The method described in Section 
(3.2.2) was initially used to determine the strain space failure envelopes for each of the ply ori- 
entations used in the laminate. Subsequently, the failure envelopes were used to develop an £2 
versus, tensile, transverse strain relationship using the progressive failure technique described 
previously. £2 was chosen as the circumferential strain Young’s moduli because as the buckling 
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Figure 3-4: Crushing Force and Cumulative Energy Absorption as a Function of Axial Dis- 
placement for the Sonar Dome. 



with respect to the laminate orientation. In reality, both compressive (bulging inward) and 
tensile (bulging outward) strains will be experienced and thus, different £2 versus, strain re- 
lationships will result. However, because the compressive and tensile strengths of the plvs are 
approximately equal (Table (3.1)). the difference is small and only the tensile £2 versus, strain 
relationships will be used to simplify the analysis. Figure (3-5) shows the £2 versus, strain 
relationship that was determined as well as the order of ply failure. 

The actual circumferential strain was then calculated as a function of x and 6 using Equa- 
tions (2.36), (2.48). and (2.44). This relationship is shown in Equation (3.9). 



numerical integration scheme to be used in Equation (2.42). 

To determine E/, a relationship between flexural strain (e/) and the axial crushing distance 



lobes are formed and bulge outward and inward, the laminate will be strained in the 2-direction 




(3.9) 



Since the circumferential strain varies both with x and <5 and £2 varies with 6. will vary 
with both x and <5. E^ was then calculated for discrete values of x and 6 in preparation for a 
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Figure 3-5: Transverse Young’s Modulus Variations with Pure Tensile. Transverse Strain for 
Frustum Section. 



(<5) was first developed using Equations (2.38), (2.39), and (2.48). (see Equation (3.10)) 




(3.10) 



Equation (3.10) was then used to determine a sj versus. 6 curve for each ply over the entire 
crushing distance for one buckling lobe (2 y/at). For each, the strain was determined at the 
midplane of the individual plys. Using these plots and the strain failure envelopes, the dis- 
placement at which each ply failed in tension could be determined. Figure (3-6) shows how this 
process was done for the 0° oriented plys in the laminate. 

Flexural stiffness (D*), E /. and laminate thickness (t) were then recalculated after each ply 
failure using Equations (3.11)-(3.13). 

[D'\ = § j\[Q)z 2 dz (3.11) 

E f = {(PT 1 ),.,} 

= nh where 



t 
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(3.12) 

(3.13) 



009 




Figure 3-6: Flexural Strain to Axial Crush Distance Relationships for the 0 Degree Oriented 
Plvs in the Frustum Laminate. 



n = number of intact plys 
h = ply thickness 

Figure (3-7) shows the variations in Ey and t as 8 increases. Notice also that as more plys fail, 
Ey becomes more erratic and diverges further and further from the in-plane modulus. This 
is expected since the laminate starts out homogenized and gradually de-homogenizes as each 
successive ply fails. 

The crushing load and cumulative energy absorption can now be determined as a function 
of axial crushing distance (Figure (3-8)) using Equations (2.42) and (3.8). Since neither Eg 
nor Ey was constant or equal during the deformation process, the simplified form of Equation 
(2.42) could not be used here. Figure (3-8) shows that five buckling lobes form as the frustum 
is crushed and that a total energy absorption of .105 MJ is realized. 
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Figure 3-7: Flexural Young's Modulus and Laminate Thickness Variation with Axial Crush 
Distance of the Frustum Section. 



0 12 
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Figure 3-8: Crushing Load and Cumulative Energy Absorption Characteristics of the Frustum 
Section. 
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3.2.4 Cylindrical Section of Nose Shell 



This outer shell of the stiffened, cylindrical section is composed a 28 ply. E-Glass (1583)/Epoxy 
(7780), woven cloth composite system laid up symmetrically 
([90.45,0, -45]-2, 90, 45, 0,0, 45,90, [90,45,0, —45)2) for a laminate thickness of 12.8 mm. Ad- 
ditional support is provided by an internal, frame structure constructed using carbon, woven 
cloth (6K 5HS satin weave IM7). composite, C-channel beams. The beams are 32 plys thick 
for a laminate thickness of 11 mm and laid-up symmetrically in the following configuration: 
[90.45.0.-45]4[90.45.0.-45]4. The shell analysis was accomplished by the same process as de- 
scribed in Section (3.2.3) and is shown in detail in Appendix (C). However, the effect of the 
stiffeners still needed to be taken into account. This was done by "smearing" the stiffeners 
into the shell and then calculating an effective thickness (t e //) using Equation (2.57). The 
detailed stiffener and t e ff calculations are included in Appendices (C) and (D) with the t e ff 
determined to be 35 mm. An effective buckling lobe length was then determined based on t e f f . 
The crushing force and cumulative energy absorption could then be determined using Equations 
(2.42) and (3.8). Figure (3-9) shows the force and energy profiles at various displacements and 
also, that the cylinder was only allowed to crush .84 m in two buckling lobes while Table (3.3) 
indicates that the section is 1.36 m long and should buckle in three lobes. This reduction in 
length was due to the effects of the transverse beam structures. The forward beam and titanium 
fittings served to force a nodal point at the beginning of the thruster door section as seen in 
Figure (3-2). The after beams and rigid titanium clevises at the aft portion of the thruster door 
create another node point and serve to terminate the deformation process. So, only the .84 m 
shell and stiffener portion between the forward and after beams wall deform during the impact. 
.49 MJ of energy are absorbed during this deformation. 

3.2.5 Analysis of Forward Hemi-head of Pressure Hull 

The fonvard pressure hull hemi-head is a metallic structure constructed from HY-80. As pre- 
viously stated, the thickness used in the analysis wall be assumed to be 19 mm due to the 
confidentiality of the actual thickness. The deformation starts as the forward anchor contacts 
the collision boundary and is pushed back into the pressure hull. The deformation is assumed 
to be of the form described in Section (2.1) and, as such. Equation (2.31) will be used to de- 
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Figure 3-9: Crushing Force and Cumulative Energy Absorption of the Cylinder Section of the 
Nose Shell. 

termine the crushing force. The flow stress was calculated from a partial stress versus, strain 
curve using Equation (3.7). The detailed calculations and the stress versus, strain profile used 
are included in Appendix (E), and the crushing force and cumulative energy absorption versus, 
displacement results are shown in Figure (3-10). Throughout this analysis, it is assumed that 
the forward hemi-head remains water tight. Further analysis must be performed to determine 
the point at which the pressure hull fails. By allowing the pressure hull to deform to the full 
axial extent allowed for the nose shell, 7.566 MJ of energy will be absorbed. 

3.2.6 Kinetic Energy Calculations 

After determining the energy absorption capability of each modeled nose fairing component, a 
determination of the kinetic energy possessed by the ASDS was needed for various speeds. The 
kinetic energy was determined using Equation (3.14) 

EE ASDS = "^{p^vehicle ^ added )^’ (3.14) 

where m ve h ic i e is the submerged mass, m a dded is the added mass, and v is the speed of the 
vehicle. 
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Figure 3-10: Crushing Force and Cumulative Energy Absorption for the Forward Pressure Hull 
Hemi-head. 
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Figure 3-11: Description of Standard Submarine Weight Terms 
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Parameter 


Weight 
(metric tons) 


NSC 


56.2 


Main Ballast Tanks 


4.8 


Free Flood 


27 


A env 


88 



Table 3.4: ASDS Weights 



The submerged mass of the vehicle was calculated using the ASDS weight report [16]. 
Weight reports for submarines are standardized and broken down into categories of weights 
based on the operating condition of the submarine. Figure (3-11) shows a brief summary of 
the standard terminology used in submarine weight analysis and how each relates to the other. 
Condition A is the dry vehicle weight. When the variable load (variable ballast tank water, 
personnel, etc.) is added, the new weight created is the normal surfaced condition (NSC). The 
main ballast tanks are filled with water to submerge the vehicle: so, the weight of this water 
added to the NSC becomes the submerged displacement (A su &). Finally, during submergence, 
non-water tight free flood areas (i.e.-between the fairings and the pressure hull) fill with water. 
This water, w'hile in pressure equilibrium with the surrounding w r ater. moves with the submarine 
as it transits, and w'hen its weight is added to A su b the resultant is the envelope displacement 
(A em; ). For this case, it is this A env that is required for the analysis. Table (3.4) shows the 
various weights used to determine the envelope displacement. 

Next, the added mass w’as determined. Since this vehicle is relatively slender, it is expected 
that the added mass in surge will not be large. However, two approximations were made using 
a sphere and an ellipsoid. The size of the sphere was set such that the frontal area was equal 
to that of the ASDS (Figure (3-12)). Once the sphere dimensions are knowm, the added mass 
can be calculated using Equation (3.15) 



thudded sphere 

(where p wateT 



1 r 2 3 
7) Pw at er^ sphere 7T ^ P water ^ sphere 



.0 metric tons 



(3.15) 



metric tons. 
1.025 3 ) 



The added mass of the vehicle v'as also approximated using an ellipsoid with the frontal area 
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19.86 m 

1.28 R m 



-2.03 m 

2.51 m 

Area = 5.17 m 



1.28 m 

9.93 m - " --- 

Side View Front View 

Figure 3-12: Dimensions Used During Added Mass Calculations 

equal to that of the ASDS (as in the case of the sphere) and a major axis length equal to that 
of the overall vehicle (Figure (3-12)). The added mass of the ellipsoid can then be determined 
using Equation (3.16) and the graphs from [17] for the added mass coefficient (C a ) 

Padded ellipsoid = C a 1 7:p water ab 2 = 2.1 metric tons (where C a = -13) (3.16) 

Since the actual vehicle shape is not as bluff as a sphere, yet not as fine as an ellipse, the added 
mass of the vehicle lies somewhere between the two calculated values. So. the average value 
was used as the vehicle added mass. Now. the kinetic energy possessed by the ASDS at various 
speeds can be determined using Equation (3.14) as shown in Figure (3-13). 

3.3 Results 

With the energy absorption of each section determined, the overall energy absorption of the nose 
structure can be determined by summing the energy absorbed by each section throughout the 
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Figure 3-13: Vehicle Kinetic Energy for Various Speeds. 

impact sequence. Figures (3-14) and (3-15) show the integrated crushing force and cumulative 
energy absorption profiles for the entire deformation process. Two curves are shown in Figure 
(3-15) depicting the energy absorption of the nose structure with and without the effects of the 
pressure hull included. Additionally, the kinetic energy curve shown in Figure (3-15) is overlaid 
on this plot. Using the curves in Figure (3-15) initial speed versus, crush distance assessments 
can be made. Without the pressure hull included, the nose structure can absorb .67 MJ of 
energy which would allow the complete absorption of all kinetic energy possessed by the vehicle 
with an initial speed of 7.5 knots. If the effects of the pressure hull are included, a total of 
8.235 MJ of energy can be absorbed. This equate to an initial vehicle velocity of 26.4 knots. 
As previously stated, the point at which the pressure hull will fail during its deformation was 
not determined: so. a the safe impact speed for the vehicle cannot yet be determined. However, 
when this point is determined. Figure (3-15) will provide this safe speed estimate. 
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Figure 3-14: Crushing Force Profile for the Entire Deformation Process. 
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Figure 3-15: Total Cumulative Energy Absorption and Vehicle Kinetic Energy. 
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Chapter 4 



Summary and Conclusions 



Throughout the course of this analysis, several new deformation models have been developed. 
The force required to deform a hemispherical shell loaded axially with a rigid plate by axisym- 
metric dishing was related to shell thickness, material flow stress (compressive and tensile), 
and axial deflection using a rigid-plastic analysis technique. Additionally, a deformation model 
for a cylindrical shell composed of elastic material, crushed axially by a rigid flat plate, was 
proposed. This model assumed a concertina type deformation mode and used shell thickness, 
cylinder radius, material stiffness, and axial displacement to determine the force required to 
crush the shell. Further, methods were set forth to include the effects of cutout sections and 
internal stiffening members on the crushing force characteristics of the shells. The effect of the 
cutout section was taken to be a reduction of the intact shell crushing force by a ratio equal to 
the shell area remaining divided by the total area of the intact shell in the local vicinity of the 
cutout. The effect of internal stiffeners was incorporated by increasing the thickness of the shell 
material to a new effective thickness determined by making the load carried by the effective 
structure equal to the combined load on the stiffeners and shell. 

These tools were used to determine the crushing force and energy absorption characteristics 
of the composite nose fairing structure and forward pressure hull of the U.S. Navy's Advanced 
Seal Delivery System in a head-on impact with a rigid wall. The complex structure of the 
ASDS was first simplified and modeled to permit the use of the tools described above. The 
shell structures were approximated as rotationallv symmetric shells (i.e.- cylinders, frustum, 
hemispherical domes) with the internal shell stiffeners being incorporated into the model in 
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two ways. The energy absorption of the longitudinal elements was captured by increasing the 
effective thickness of the shell in the vicinity of the longitudinal beams, while the transverse 
stiffening elements served only to force nodal points in the shell buckling pattern since their 
energy absorption during a head-on collision is minimal. Further, major internal components, 
such as the titanium clevises (used to connect the internal beams) and the forward anchor, were 
treated as rigid, massless objects. 

The structure was divided into a hemispherical sonar dome, a frustum and stiffened cylin- 
drical nose shell, and a hemispherical forward pressure hull based on the shape, material com- 
position. and thickness of the actual structure. Using laminated plate theory and a progressive 
failure model, stiffness versus crush distance and material thickness versus crush distance rela- 
tionships were developed for each of the material systems used. These relationships were then 
used in conjunction with the rigid-plastic hemispherical and elastic cylindrical shell models 
to predict the force deflection and energy absorption of each component. The component- 
level energy absorptions were then summed to approximate these characteristics for the overall 
structure. Finally, the absorbed energy was equated to the vehicle’s pre-impact kinetic energy 
in order to determine an initial speed just prior to the collision. The composite nose fairing 
structure was deemed capable of absorbing the impact energy associated with a vehicle speed 
of 7.5 knots. If the effects of the pressure hull were included, the structure could absorb the 
kinetic energy from a 26.4 knot impact. 

4.1 Areas Additional for Study 

This analysis has taken the first step in determining the energy absorption capability of the 
ASDS. However, considerable work remains to be done before a final assessment can be made 
of the vehicle's overall crashworthiness. The following is a list of topics for further study in this 
project: 

• Model validation and calibration: The models developed herein have not been cali- 
brated or confirmed via actual test data. Cylindrical and hemispherical composite sections 
of the appropriate length to radius and thickness to radius ratios need to be constructed 
for each composite material and crushed under quasi-static and dynamic conditions. The 
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results should be used to confirm, improve, and calibrate the methods developed in this 
thesis. 

• Rupture Point of Forward Pressure Hull: The overall objective of the project is to 
determine a safe impact speed for the vehicle such that during an impact the pressure hull 
will not rupture and the crew will survive the collision. The rupture or collapse point of 
the forward pressure hull needs to be determined and may vary with depth of the vehicle. 

• Strength Analysis of Nose Fairing Mounting Plates: The composite nose fairing 
is mounted to reinforced metal plates that are welded to the forward hemi-head of the 
pressure hull. An analysis needs to be done to determine whether these structures will be / 
capable of withstanding the force of the impact or if they will bend or shear off completely. 

• Energy Absorption of Internal Components: The energy absorption potential of the 
internal components (such as clevises, angle bulkheads, the anchor, the sonar equipment, 
etc.) have not been included in this analysis. The effects of each of these components 
can be determined individually and added to the curves developed in this report to get a 
more accurate approximation of the energy absorption of the vehicle. 

• Effect of Impact Angle on Energy Absorption: This analysis looked at an "impact 
sequence" for a head-on impact angle. In reality, this situation is likely never to happen. 

A relationship should be developed for various impact angles using model testing or other 
means to relate the reduction in energy absorption that is realized by an off-axis impact 
to that of an orthogonal impact. 
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Appendix A 



Sonar Dome Analysis 
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Ply Stiffness and Layup Inputs: 



1 . Material System: Kevlar (K285)/ epoxy (7714) woven cloth 

2. Layup Configuration: (90, 0 ) 5 , 0,0, (90, 0 ] 5 



n 0 :=12 


0 

11 

z 


n 90 :=io 


N p45 : =0 



3 Total Number of Plys Totai_piys:= n 0 +n n45 +N p^+N 90 

4. Ply Thickness h :=.oi in h mbber := 15 in 

5. Laminate Thickness t L :=Tota]_piysh-i-h rubber 

6 . Strain Applied to the Materia! (circumferential): 

P=atan(s y /£ x ) p :=o-deg 

Shear Strain: e s :=o 

7. Number and Orientation of Ply Angles: 

Plys: =4 
a : = 1 .. Plys 
c ; = 1 .. 3 Plys 

8 Engineering Constants: 



Kevlar 




Rubber 


E X :=4 7-MS1 


E s : = 45 MSI 


^ rubber • 1 304 -psi 


E v :=4.7-MSl 


E m :=3 4-GPa 


v rubber >- 499 


v f :=5 


G m :=1.26 GPa 




v m •” 5 


v x :=]5 





■deg 



t L = 9.398°mm 



/ 
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Model Parameters 



1. Radius of Base of Sonar Dome. 

2. Axial Depth: 

3. Radius of Curvature 

4. Radius to Center of Cutouts 

5. Axial Distance to Center of Cutout 

6. Slope of Dome Surface at Cutout Center 

7. Diameter of Cutout 

8. Radius where Cutout Starts 

9. Radius where Cutout Ends 
10 Number of Cutouts 



R model <_ 36 *n 



D model *"‘5 in 



. ,_1 D model + R model 

1 curvature '~—'~ 



el 2 ' 1 
D model 



re— */i 



cutout • K curvature - >\7 K curvature — K cutout u cutout = 



-D cutout curvature 



cutout cutoufh^ R curvature 



) =58.003°deg 



d cutout •“* 1 ' n 
r 0 :=23.693in 
r f ;=29.66 m 

n cutouts ,- 3 



Ply Strength Parameter Inputs: 

1. Tensile Strength in Longitudinal Direction 

2. Compressive Strength in Longitudinal Direction 

3. Tensile Strength in Transverse Direction 

4. Compressive Strength in the Transverse Direction 

5. Shear Strength 



X:=56KS1 



X prime' 14 R S1 

Y :=56 KS1 

Y prime ,-14 R S1 



S :=4 5-K.SI 



Matrix and Fiber Damage Inputs 

E f_star :=01 t|y :=5161 n!=1 

E \v star 1 " 



- m_star 



:= 15 



n s := 316 
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Damage Ply Calculations: 



p - Ex r - 

E fy — Cj (V • 



E y_damage<f E ) • — 



E s_damaged* E * ’ 



i+-n s 



Building Ply Data: 

Undamaged ply, ply with m 



, v m 


v m 1 
" s — 


J 

E fc := 


. v m v m 

+n y~ — 




E s 


G m 


ty 


E y E m 




] + 1 


v m 

n >— 


-j-1 




v m 


E fy F 


E m 












- 








v m 


-1 






1 1 


nS ~ 






v m 




G m 







Undamaged ply, ply with matrix damage and ply with fiber damage. 



‘ E x ‘ 




E v 




E S 


E x 


E y y : = 


E y_damaged E m_star 


E ss 


E s_damaged E m_star 


E x E f_star 




E \_damaged E f_star 




E s_damaged E f_star 





v X 




E xy_star 


v XX 1 


v x E m_star 


EE xy_star' 


E xy_star E m_star 




v x E f_star 




E xy_star'E f_star 



xx nrim£U :=x 



primeq ,-A prime 



Rotation Matrices: 



-yy q 



^'a 

4 

n a 

, 2 '"a 2 

i 2 n a 2 

'a\ 

„ 3 



n 4 


2- ml 


n 


2 


4 ml 


a 


a 


a 




a 


ml 4 


2 ml 2 - 


n 


2 


4 ml 


a 


a 


a 




a 


ml a 2 "a 2 


m, a V 


n a 


4 


-4 ml a 










r ”> 


ml a ' "a ‘ 


-2 ml a - 


n a 




[""a 



m, a 

ml 



n ml ^ n 2 [ml • n J — 
a a a |_ a a 

( 



ml a 

2-1 ml J n — ml • 
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Ply Stiffness Calculation: 

Calculates the zero direction stiffness matrix for the various ply states {undamaged, matrix damaged, and 
fiber damaged) 



E « q 'yy q E x.Vq n 


^ rubber v rubber ^ rubber n 


2 2 u 

1 — v rubber ' ~ v rubber 

v rubber E rubber E rubber n 


1 “ v x\j v y> q 1 “ v xxq v >7q 
'\.\q ^y\ q ^ X>q ^ Q rubber - 


2 2 
' — v rubber ' " v rubber 

n n 2 E n,bber 


1 ~ v xxq v yy q 1 - v x\q v yy q 

0 0 E,, 

L qj 


1 i~ v rubber 



Creates the stiffness matrices for every combination of ply angle and ply state. Qnew is a 3 x 4 matrix 



Q rol q , a 



: -Rotation l a 



r On 1 






q 1,1 






Q q 2,2 




^ rot q,a | ^ rot q,a 3 ^ rot (],a 5 


°q 1,2 


0 new q a : 


Q r °t q ,a 3 Qr °tq,a 2 ° rot q,a 6 


°q 3 , 3 _ 




Q rot q , a 5 ° rot q,a 6 Qrot q,a 4 



Stress Space Strength Parameters: 

A separate F must be calculated for each ply state 



s 



i 



xx q xx pnmeq 

xy_star q 




0 



^XX q .XX prjmeq .YY q YY pnmeq 

1 



0 



0 



0 




F 



1 1 

xx— 



YY~“ YY 

q 



j 

prime q 
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Strain Space Strength Parameters: 

1 . Strain strength parameter for each ply state 



Su[V,] : ' 


F \ 1,1 1,1 Q q 1 ,2 


0 


Fa q i,: Q q i.i Q q 1,2 ' 


+ F *q 1 ,2 [ Q( i 1 , T Q q 2,2^1 Q q 1,2] ] 






+ F ^ 2,2 V2 Q < 2,2 




S i,i' Q q i,i' Q q 1,2" 


Fa q u [ Q q 1,2] 


0 


F *q i ,2 [ Q q i.i °q 2.2 + [ Q q 1,2] ] 


+2 F ^ 1,2 Q( i 1,2' Q q 2,2 




+ °qi.2 Q q:.2 


+ S 2 . 2 [V 2 ] 2 




0 


0 





F b q ,0, , F b q , Q q , . 

% 1 Q q I.; + F »q ; °q;.: 



2 Calculating G matrices for all combinations of ply state and angle: 





G base! , 

G basel | ^ . 




Rotation 1 ( 


G base 1 1 j 2 
G basel j _ „ 
G basel 2 j j 
G basel 2 




Rotation 1 , 






1 


G basel 2 | 1 

G basel t , 

- j,3 

G baselj | |" 
G basel- 




Rotation 1. 






1 


G basel 3 j 0 

G basel 3 , 

->.3 

Cbasel, 

G basel] 





G 2irans :: 
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Rotation^ 



RotationZj 



G base2 t ' 
G base2j 
G base22 
G base2 2 2 
G base2 3 
G base2 3 
G base2i 
G base2 j 
G base2 2 
G base2 2 
G base2 3 
G base2 3 
G base2j 
G base2j 
G base2 2 
G base2 2 
G base2 3 



G 1 trans 1 ' 



1 "> ^ 



Rotation l 0 


G basel i j t 




G base 1 , 3 3 




G base 1 2 j j 


Rotation 1^ 


G basel 2 2 

G basel 2 j t 

G basel -> , •> 

“ j,3J 

G basel 3 u 


Rotation^ 


G base 1 3 :>2 
G basel 3 j t 
G base 1 3 3 3 
G basel i j j 


Rotation Ij- 


G basel j t 
G base^ ] ^ 




G basel i 3 3 
' G basel 2 , ,1 


Rotation U 


G baseH 0 -> 
G basel 2 j t 




G baseH 3 3 
' Gbase^ K1 1 


Rotation ly 


G basel 3 t t 
G basel 3 12 
G base 1333 

’ G basel 1 j f , 


Rotation^ 


G basel 1 t t 
G basel j t 
G base li , , 

3 » 3 J 

G basel2 j j 


Rotation^ 


G base 1 2 t 2 
G basel 2 12 




G base >233 
G basel 3 , , 


Rotation!^ 


G basel 3 t 2 
G base 1 3 j t 




G basel 3 3 3 



Rotation^ 



Rotation2 4 



Rotation2 4 



Rotation2_j 



3 1 

G base23 2 
c base2, ^ 
G base2] ? 
G base22 j 
G base2 2 2 
G base2 3 ] 
G base2 3 ? 
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■ G 


Itransj. ( 


G Itransc 3 


G ltransj. . 


G 2transc 


G C : = 


G 


Itransj. , 


G 1 transc 


G Itran^ 6 


G 2trans c 




G 


Itransj- . 


G ltransc 6 


G ltransj. 4 


G 2trans c 



Formation of Strain Space Failure Envelopes: 

1. The following relationship is used to determine the failure envelope: Gj -EjCj-i-Gj**! 

2. After shifting to polar coordinates, the above equation can be solved for r using the quadratic equation: 

E| = r-cos<) E-, = r sin<t> A f -|-B r— 1=0 

3. Calculation of coefficients for the quadratic formula 

A(4».b) := G b ^ ^sin( <>)"-(- 2- G b ^ ^ -cos( <t> ) sin( <{>)-!- G b ^ ^ cos( <|> )*" 

B <t> ,b,£ s :=2- G b ^ , sin(((») E s -(- G b ^ ^*sin< <J> G b ( ^ ^os( <> ) -e s -(- G b ^ ^ cost <b ) 

C <t> ,b,E s : = G b 3 4 -e s -h G b 3 3 -e s “- 1 



-B <>,b,E s +JB ()),b.E s " — 4A((t>,b)C <|>,b,E s 

r«ii,b.E s i- 1 — 

5 2 A( <> ,b) 

Strain Space Failure Envelopes: 

4> :=0deg, I deg.. 360-deg y(x) :=tan(P) x x : =0,1 ^ 1— .. cos( P) e x <|> , b, e s :=r <t>,b,E s cos( <t> > Ey6,b,E s i=r t,b,E s sin( <}> ) 



Ply Strain Strength Envelopes 




0 ply -Undamaged 

0 plv-matrix damaged 

90 ply-Undamaged 

90 plv-matrix damaged 

45 ply-Undamaged 
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Calculation of Degraded Laminate Stiffness: 

The first ply failure strain is first calculated. Subsequently, a progressive failure of plys is allowed and the 
stiffness calculated based on the ply damage history. 



plotl - = 



al<-l 
bl«- 1 
count*- 2 
nl*— 1 
s*-l 
while 1 
cl | *— 0 

CTl,«-Ot 



z x p,al ,z s 
z x P,bl -|- 3 , z s 
t 'tan( p ) 

1 



C0Um M M h ^ bber 

N 0 S + N90 S-I ^ 

count 4- £[ A star count J i , ij 

’z\ 



count 

count 



al«-al -j-2 if e! =e x P.al,c s 



■ Q- Q ne\v a | ( | 'S-f-N 99- Q new^i 2 s_ h fj ^ rubber j 



bl*-bl -t-l ifsl coun = 
count*— count-)- 1 
break if a 1 b 1=6 
break if al >3 
break if bl >3 

ct! 



: x P»hl +3 ,e s 



tempstore*- augment el 



MPa 



• 4 236 1 

z rubber “ | 






" 


z rubber) ,j 


iii ! — 1 rows z ^5^ 


CT rubber, jj 


Q rubber 


0 










0 J. 



plot! -stack plotl, augment z m bber- 



CT rubber 
MPa 



Tensile Flow Stress Determination curve* x) :-linterp plot 11 > ,plot <2> ,x 



o 



r pio^ 1 



. 1 > 

P |ot rows* plot) 



rows* plot) 



curve* x)dx MPa 



ct , = 53 . 354 °MPa 
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Calculation of Degraded Laminate Stiffness (Compressive): 

The first ply failure strain is first calculated. Subsequently, a progressive failure of plys is allowed and the 
stiffness calculated based on the ply damage history. 



plot 1 .'= 



p ! = 180deg 
al«— 1 
b1<- 1 
count«-2 
nl < — 1 
s< — I 
while 1 
£^<-0 
ol j«-0-MPa 

i < — a 1 



x P * a 1 , e j 
P * bl -j- 3 , e s 
tan( p ) 



1 star count 



k, vi ^ rubber 

N 0 S + N 90 u 



[no 



/"v v i ~ h rubber 

' Qnew a j j s TN9o- ^new^j ^ s + j; Q rubber 



^ 'count 4- |J A star count J],]J 



v star count 



count 

^count 



al«— al -j-2 if El coun =E x P,al,E s 

''count - 
counu— count -j- 1 
break if al-j-bl^ 
break if a 1> 3 
break if b 1 > 3 

rr - 1 

tempstore*- augment e! 



MPa 






iii :~1 .. rows e m5ber 



a rubber^ *” 


Q rubber 


E rubber,;; 
0 






0 J_ 



plot I— stack plot I, augment e rubber- 



° rubber 
MPa 



Compressive Flow Stress Determination curvet x) l=linterp plof* > ,plot < “ > ,.\ 



plot 



ro\vs( plot) 



rows( plot) 



curve( x)d x-MPa 



a c = 8 41 °MPa 
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Force vs Axial Crush Distance Determination 



nn :=30 bbbl = l..nn 



1. Establishing Displacement Array _ ._ D model -(bbb-i) 

6 bbb- ; 

2. Calculate Open Section Correction Factor 

a( 5 ) 2 R cur vature~ 6 





[ r o + d cutout c 0 * 10 )- a* 6 ) A r 0 -t-d cutout cos( o ) - a( 6 ) 


cost u> ) tul " ul cost 0 ) ) 


a<5) 



n cutouts^ 5 ) 
e(o) .-1 

2 71 

3. Calculate Unbalance Flow Stress Correction Factor 

a := T CT --. a = 1 728 

a c +a , 

4. Calculate Crushing Force for Displacement Array 

1 3 

2 ' l L 2 ' 6 bbb 

Pl(xx) :=(lmterp(5.P,xx)) 

Energy Absorption of the Structure: 

Energy( 51 1 ) :=linterp( 5 ,P,61 1 ) 

P 6 bbb 

Cumulative_Energ\^ bb ! = J Energy( x)d x Cumulative_Energy row ^ g ^ = 0.058°MJ 





1 3 1 


■ 


1 3 r 


p bbb :=lf 


. a6 bbb ^o* 4 -"V2-o c -o 2 -t L 2 - S bbb 2 .«f 


_ a 5 bbb 5 bbb ' 


r - - ~ 

4 -7:^/2 <r c a“ t L • 5 bbb “ 
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Appendix B 



Frustum Section Analysis 
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Ply Stiffness and Layup Inputs: 



1. Material System: E-Glass (1583)/Epoxy (7780) woven cloth 

2. Layup Configuration: [SO.AS.O^Sk.OO.AS.O.O^S.QO.fOO.AS.O.^Sfe 



n 0 :=6 N n45 :=4 

N 9 o : = 6 Np4 5 :=6 



3. Total Number of Plys Totai_pivs:= n o+n n45 +N p 45 -hN 90 

4. Ply Thickness h := oi s in 

5. Laminate Thickness ti :=Totai_pivsh 

L t L = 0.010l*m 

6. Strain Applied to the Material (circumferential): 

p=atan(£ y /E x ) p :=89.5<ieg 

Shear Strain: e s :=o 

7. Number and Orientation of Ply Angles: 

Plys! =4 
a :=l .. Plys 
c : = 1 .. 3 Plvs 



■deg 



8. Engineering Constants: 



e x :=3 MSI 



e s :=.55MS1 



E y : = 2 7 MSI 



E m ,= 3 4-GPa 



G m : = I.26GPa 



90 

45 

0 



-45 

90 

45 

0 

-45 



sequence \ - 



90 

45 

0 

0 

45 

90 



-45 



0 

45 

90 



-45 



0 

45 

90 



v f :=5 

v m := 5 

v x := 15 
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Component Parameters: 

1. Radius of Shell: 

2. Length of Shell: 

3. Wave Length of Elastic Buckling Wave 

4. Number of Buckling Waves Along Shell 



^ shell ' 3 -25 ft 
L shell : =31 -ft 



h:=^ 



Number waves 



._E shell 

~ m - 



Ply Strength Parameter Inputs: 

1. Tensile Strength in Longitudinal Direction 

2. Compressive Strength in Longitudinal Direction 

3. Tensile Strength in Transverse Direction 

4. Compressive Strength in the Transverse Direction 

5. Shear Strength 



Number wave5 . = 4 733 



„ E shell 

efT ‘ 2 Number wa . 



X:=39 1KS1 
X prime ’ = 33 8-KS1 
y:=39.iksi 



Matrix and Fiber Damage Inputs 

E f_star :=01 n y := 5161 n := 1 

p -i- „ •-ha F w star - - 5 

E m star • - 1 5 n s >_316 

Damage Ply Calculations: 






1 n s — 



E fv :: 



1 -»-n V — 



E \ _damagecf E ) • - 



>_damagecf E ) • - 



1 


i l 


v m 1 

n y 

y Vf 


, v m 

1 -f-n y — 


[E fy + F 


Em J 






v m 


1 


1 1 


nS “ 


. v m 

i -f-n s — 
V f 




G m j 
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Building Ply Data: 

Undamaged ply, ply with matrix damage and ply with fiber damage 



x c m star 



xy star & m_star 
^ xy star E f star 





' E x ■ 




Ey 




II 

Si 

U1 


Ex 


E yv • 


E y_damaged E m_star 


E ss •“ 




_E x E f_star_ 




E y_damaged E f_star 






v X 




E \y_star 


XX £ 



E s_damaged E m_star 
E s_damaged E f_star 



' primeq * ^ prime 



YYq •-'t YY p nmec| • Y pnme 



' y>fl XXn 



-xx q 



Rotation Matrices: 



ml a 4 


"a 4 


2 • ml a 2 - 


n a 


2 


4 


"'a 2 V 


"a 4 


ml a 4 


“ m *a 2 


n a 


2 


4 


m 'a 2 n a ~ 


"'a 2 "a 2 


m 'a 2 n a 2 


m 1 a 4 -i" 


n a 


4 


-4 


"'a 2 "a 2 


-a 2 "a 2 


m *a 2 n a 


-2 ml a 2 




2 




ml a 3 "a 


- ml a "a 3 


ml a "a 3 - 


m, a 




[ ml a* 


"a 3 " ml a 3 "a] 


m, a "a 3 


- ml a 3 "a 


m *a ^ n a — m 'a 


"a 3 H 


[ ml a 


3n a- ml a "a 3 ] 
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Ply Stiffness Calculation: 

Calculates the zero direction stiffness matrix for the various ply states (undamaged, matrix damaged, and 
fiber damaged) 



E XXq 


v yy q ' E xx q 


>- v xx q v » q 


| - v .\x q ' v >y q 


vx v E »q 




, — v *Xq' v »'q 


l-' Mq v v; q 


0 


0 



Creates the stiffness matrices for every combination of ply angle and ply state. Qnew is a 3 x 4 matrix 



Qr0t q,a :=R0tatl0nl a 



Q q 

Q q 

Q q 



1,1 

1,2 

3,3 



rot q,a , 


QroI q.a3 


^ r0I q,a 5 


rot q,a 3 


Q ”»q,a : 


QroI q.a 6 


rot q,a 5 


Q,, ^H 


S.»4 



Flexural Stiffness: 

1. Indexing variables ff :=i ..Total _piys-j-i 

eg :=1 ..Total_plvs 

2. Establishing an array of thickness positions for each ply 



z fT ' 



:=llL-Hfr-l).h 



3. Assigns the appropriate stiffness matrix based on the layup sequence 



Qflex 



gg 



Q new, j if sequence gg =0-deg 
Q new, ^ if sequence gg =90-deg 
Qnew, 3 > f sequence gg =45deg 
Qnew, 4 otherwise 
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Stress Space Strength Parameters: 

A separate F must be calculated for each ^ly state 



' x y_star q 



XX 1 xx pnme q ^xx q XX prlm5q Y^' q YY prllne(] 



/.\^pr, m5q VY q VY pnmC() 



' 1 1 

" VA q ' VA pnmeq 
1 1 

11 'q 11 pnme q 



Strain Space Strength Parameters: 

1 . Strain strength parameter for each ply state 





Su[Vi] 2 ' 


Fa q 1,1 Q q 1,1 Q q 1,2 




0 




*' Fa q 1,; Q q l,r Q q 1,2" 


+ F Aq 1 ,2 [ 0( 1 1,1* °q 2,2 + [ 


Q < I.:] 2 ] 






+ F =q 2 , 2 [ Q q,. 2 ] : 


F ^ 2,2 1 , 2 ' Q q 2,2 






G basel q 


Fa q i.i Q q 1,1 Q q 1,2 


F ®q 1,1 [ Q q 1,2] 




0 




F *q 1 ,i[ Q q 1,1 °q 2.2 + [ Q q 1,2] ] 


+ 2 F a 0 O 

^ 1 2 1,2 v q 


2,2 






F a O Q 

^ 2,2 ^ 1,2 q 2.2 


+ S 2 . 2 '[V . 2 ] 2 








0 


0 




?a * 3,3 [ 3,3] 



F b q i Q q i,i + F b q : Q q 1,2 
F b q , Q q |. 2 + F »q 2 '°q 2 ,2 
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2. Calculating G matrices for all combinations of ply state and angle: 



Rotation 1, 



Rotation lj 



Rotation 1 1 



Rotation 1-, 



Rotation!-, 



Rotation 1-, 



G 1 trans : = 



Rotation 1, 







G base 2 ] 




Rotation^ 








G base 2 , 






G base 22 l 




Rotation2| 




Gbasel, , , 








G base 2 2 


G basel] 








G base 23 








Rotation2| 




G base 1 1 j 








G base 23 


G base 1 1 3 3 








G base 2 ] . 








Rotation2-, 




° base| : i.i' 








G base 2 ] 


G basel 2 








G base 2 2 . 








Rotationij- 




G base 1 2 j 








G base 2 2 , 


G baselT , , 








G base 23 


’ 






Rotation 22 ‘ 




G base 1 3 , 








G base 23 2 






G 2 trans ' 






G base 1 3 - _ 








' G base 2 , 


G basel 3 ^ 






Rotation2 3 • 


G base 2 ] 










G base 2 2 . 








Rotation2j 




Gbasel, , 








G base 2 2 


G base 1 j , 








G base 23 


^ ’ * 






Rotation 22 • 




G basel| j ^ 








G base 23 ^ 


G basel , 33 








G base 2 ] 








Rotation2^ • 




G basel 2 | j 








G base 2 ] 


G base 1 2 








' G base 2 2 








Rotation2^- 




G basel 2 j 








G base 2 2 ^ 


G basel 2 3 3 








G base 23 j 


’ - 






Rotation2^ 




G basel 3 ( j 








G base 2 3 , 


G basel 3 








G base 1 3 ^ ^ 








G basel 3 , , 








G basel j ' 








G basel j 








G baselj j 








G baseli „ , 
1 
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' 3 . 3j 





L 




' G base 1 2 , ,] 




G baseH t 


Rotation^ 






G basel 2 j 




G basel 2 3 3 
G basel 3 u 
G basel3 ^ 


Rotation 1 ^- 




G base 1 3 , , 




G base 1 3 3 3 




G basel i j ( 




G basel i 


Rotation^ 






G base 1 1 j , 




G base 1 ! 33 




r G basel*. , ,1 




-1,1 




G basel 2 n . 


Rotation!,- 




4 


c basel 2 , ^ 

G basel-) „ , 

- j, 3 J 

‘ G basel 3 
G basel 3 


Rotation 1^- 






° basel 3 1.2 

G basel3 3 3 



ltransj. j G ItranSj. 3 


G ltrans c 5 


G 2trans c 


ltransc , G ltransj. 


G ltran Sc 6 


G Ztransc 


ltransc . G ltransj. ^ 


G ltransc 4 


G 2trans c 
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Formation of Strain Space Failure Envelopes: 

1. The following relationship is used to determine the failure envelope Gj ■ z t ej+Gj e=i 

2. After shifting to polar coordinates, the above equation can be solved for r using the quadratic equation: 

Ej = r-coso e n = r-sin<t> Ar fBr-R 

3. Calculation of coefficients for the quadratic formula 

A(4>,b).‘= G b ^ ^sin(«t»r-|-2- G fe ^ ^cosi $ ) sin( <j> )+- G b ^ ^ cost <> ) 2 

B <J>,b,E s :=2 G b 2 3 sin(<|>) e s -|- G b o 4 sin( 4> )-|-2 G b ( ^ cos( <J> ) e s -|- G b j 4 <os( ) 

C <t>,b,E c ■“ Gl. •£ c+ Gv -e c — 1 
s b 3.4 b 3.3 5 

_-B b.b.E s -)-jB <j>,b,£ s “-4A($,b)C 4>,b,E s 

r o ,b,E s :=. z 

s 2 • A( i)> , b ) 

Strain Space Failure Envelopes: 

4> :=0-deg,l deg..360deg v( x) :=tan{ (J )x x 1=0. .. cos( (3) E x $,b,e s t=r <j> ,b,E s cos( 4> ) e y 4> , b , e s :=r<>,b,£ s sin{ 4> ) 



Ply Strain Strength Envelopes 



i nmi ! 






\ /A 


, K. . 


-do4 -d 03 -J.02 j 4djoi 


0 O'OI // 0:02 '* N 0'03 

jy' x 











Longitudinal Strain 

0 ply-Undamaged 

0 ply-matrix damaged 

90 ply-Undamaged 

90 ply-matrix damaged 

45 ply-Undamaged 
45 ply -matrix damaged 

45 piy-Undamaged 

-45 ply -matrix damaged 

Strain Relationship Line 



78 



Calculation of Degraded Laminate Stiffness: 

The first ply failure strain is first calculated. Subsequently, a progressive failure of plys is allowed and the 
stiffness calculated based on the ply damage history. 



plot := 



a 1 < — 1 
bl«- 1 
cl«-l 
d 1 < — 1 
count*- 2 
n 1 * — 1 
while 1 
cl 



ol | <— 0 MPa 

c x (J,al ,e s 
e x p.bl -)-3 .e s 
e x P,c1-(-6.e s 
e x p , d 1 -|-9 ,e s 

"count 4-8 ' count tan ^ ^ 

A star. 



count 4- 



r count 4— N 0 tN 90 tN n 45 tN^ [ N °' Qnew al,l + N 90' Q new bl >2 n45‘ Q new cl 3 + N p45' 

^ y ' count 4- A star count ] 2 2 j 







el count 


a ' count 4- 


A star count 


f2 

count 
e s 



al*— al -j- 1 if el count =e x P,al,E s 

bl*-bl + 2 if el count =€ x P.bl-h3,e s 

Cl4-cl-|-l if El coun =c x P.cl-h6,E s 

dl*-dl + l if El coun =e x P.dl + 9,E S 

count*— count-j- 1 

break if al +bl -f-cl -(-dl=12 

break if al>3 

break if bl>3 

break if c 1> 3 

break if dl>3 



tempstore*- augment e2.— _ 
E vl 

augment tempstore.^- 



nevv dl ,4 
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Ply Transverse Young's Moduli as a function of Strain 



Circumferential Ply Stiffness (Tensile) 




0 


0 


0 


0014 


247 801 


16 449 


0015 


180.252 


10.741 


0.016 


126.253 


7.07 


0.025 


85.233 


2.765 


0.171 


94 954 


0 554 



Transverse Stiffness Variations with Axial Location and Axial Crush 
Distance: 



1 . Setting up stations along the x* axis and axial crush distances at which E y is to be calculated 



= 1.17 zzz :=l..n 



= _--L- H efr 



MO' 



v-3 



r3 



9 10' 
0.011 
0.014 
0 02 
0.029 
0 041 
0 068 
0.093 
0.1 
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2. Calculation of circumferential strain based on axial position and crush distance 



._ x pnme^-o tempyvy 
0 zz2, R shell 14 etT 



I HefT ) 

^ tempyyy ^ 



3 Assignment of the appropriate Young's moduli to the correct axial position and crush distance based on the 
transverse Young's Moduli vs strain plot and circumferential strain. 

0 GPa if e >171 

°zzz, yyy 



2 765-GPa if 025>e o ^016 

0 zzz,yyy 

7 07 GPa if 016 >e a >015 

w zzz,yyy 



16 449 GPa otherwise 



plot = 



0 


0 


0 


0.014 


247 801 


16 449 


0 015 


1 80.252 


10 741 


0.016 


126.253 


7.07 


0 025 


85.233 


2.765 


0.171 


94 954 


0.554 



Flexural Strain and Stiffness Calculations, and Thickness Variations: 

1. Flexural Strain as a function of crush distance oi8m = 4.572*io' 4 «m 

_£ 5:=0..001..2^R shell t L 



E xf( 6,gg):=2-^| +I gg 






Flexural Strain*0 Orientation Ptvs 




PIv 5 
- Ply 9 

• Tensile Failure Strain for 0 Oriented Ply 
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l-lexural Strain Flexural Strain 



Flexural Strain-45/-45 Orientation Plys 




Plv 4 

Plv 6 

Plv 8 

Ply 10 

Tensile Failure Strain for 4S/-45 Oriented Ply 

--- Compressive Failure Strain for 45/-45 Oriented Ply 



Flexural Strain-90 Orientation Plys 




Plv 7 
Ply 1 1 

Tensile Failure Strain for 90 Oriented Ply 

- - - Compressive Failure Strain for 90 Oriented Plv 
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2. Determination of flexural stiffness and thickness variations with crush distance This a manual operation where 
the plys are failed progressively as crush distance is increased to 2H using the graphs above and the included ply 
strain limits. The thickness is determined by removing ply thicknesses as they fail. 




D star 



12 1 

’ (Tio 13 



Total_plys 

Z Q nex, x ([Wi)] 3 - 

xx= 1 



E xf • dstarjj 



E xf = 20 6843°GPa 



0 



.002 




‘164811* 




'22 h' 


003 




164811 




22 h 


.005 




17.3293 




19 h 


.006 




15.9983 




15 h 


.007 




16.4377 




13 h 


.008 




18.3322 




1 2-h 


009 




16 5229 




1 1 -h 


Oil 


H efif 


16.7947 




10 h 


.014 


H ^ flexure " 


17 0984 


■GPa t 


9-h 


02 




14 2239 




8h 


.029 




16 142 




7-h 


041 




17.2425 




6-h 


.068 




18.3592 




5-h 


093 




12 799 




4 h 


H 




18 4517 




3h 


TrT 




20 6843 




2-h 


2-H 




20 6843 




2-h 



m 
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Force vs Axial Crush Distance Determination: 



t L 2 H eff E flexureyy^,'^ shell' 1 fyyy 

2 2- E0 ppp,yyy X P nme ppp n ^ ~ 3 



R shell H eff ppp= 1 



Energy Absorption of the Structure: 



Energy! 611 ) .'=linterp< 5.P.611) 



Cumulative_Energ>^^, 




Energy! x)dx 



Cumulative_Energ> row ^ g ^ = 0 02 1 °MJ 
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Appendix C 



Cylindrical Section Analysis 
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Ply Stiffness and Lay-up Inputs: 

1 . Material System: E-Glass (1583)/Epoxy (7780) woven cloth 

2. Lay-up Configuration: [90,45,0,-45]3,90 I 45 I 45 I 90 I [90,45,0,-45]3 



n 0 :=6 

n 90 :=8 



N 



p45 1 



3. Total Number of Plys 

4. Ply Thickness 

5. Laminate Thickness 



Total_ply S :- N o-1-N n45 -|-N P 45-J-N 90 
h : = 0I8 in 
t l : = Total_plys h 



6. Strain Applied to the Material (circumferential): 

P=atan(e y /£ x ) p :=89.5-deg 

Shear Strain: e s :=o 

7. Number and Orientation of Ply Angles: 



Plys: =4 
a :=I .. Plys 
c : = I .. 3-Plys 



■deg 



8. Engineering Constants: 



E s :=.55 MSI 
E m :=3 4 GPa 



G _ : = 1.26 GPa 



0 

-45 

90 

45 

0 

-45 

90 

45 

0 

-45 

90 

45 

45 

90 

-45 

0 

45 

90 

-45 

0 

45 

90 

-45 

0 

45 

90 
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Component Parameters: 

1. Radius of Shell: 

2. Length of Shell: 

3. Wave Length of Elastic Buckling Wave: 

(shell only, no stiffener effect) 

4. Stiffener Parameters 

Undamaged Stiffness 

Number 

Area of Single Stiffener 

5. Open Section Parameters 

Reduction Parameter 



R shell -“3.5 -ft 
L shell : =33 m 

H shell 1 L 

E carbon 5j.793-GPa 
N stiffener’ - ® 

A stiffener := ( 5 52 6 -in-|-( 8 in- 2 45in))- 45in 
S :=62 5-deg 

Y y=0 653 



Ply Strength Parameter Inputs: 

1. Tensile Strength in Longitudinal Direction 

2. Compressive Strength in Longitudinal Direction 

3. Tensile Strength in Transverse Direction 

4. Compressive Strength in the Transverse Direction 

5. Shear Strength 



X:=39 1 -K.SI 
X ^ime : =33 8-KS1 



Y:=39.1 KS1 
Y prime 8 KS1 
S :=9.41KS1 



Matrix and Fiber Damage Inputs 



E f_star 01 


ti v :=5161 


n 1 


^ m_star ,_ 1^ 


t \ s := 316 


F xy_star 
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Damage Ply Calculations: 



E fx G fx :: 



1 v m 
i-m s — 


v m 

1 S" — 


-1 


. % m 

1 +•’1 v 


v f 


v f 


E a, 


' v f 


E s 


G m 


E fy 


L E y 



E y damaged* E) •“ 



E s_damaged* E ^ • 



i 


1 + 


v m T 

71 >' v f 


V m 


E fy F 


E m 


i -l-n y — 








v m 

T) s 


l 


1 1 


Vf 


, v m 

’-M s — 


G fx^E 


G m J 



Building Ply Data: 

Undamaged ply, ply withTmatrix damage and ply with fiber damage. 



‘ E x ' 




E > 




E s 


E x 


ll 

LU 


E v_damaged E m_slar 


Ess : = 


E s_damaged E m _star 


E x E f_star 




_ E y^damaged E f_star 




E s_damaged E f_star 





V X 1 




E xy_star 


v xx •“ 


v x E m star 
v x E f_siar 


EE xv_star“ 


E xy_star'E m_star 
E xy_star E f_star 



x primeq 1 prime' 



YY q :=v YY pnmeq :=Y pnme 



>V xx q eT 



Rotation Matrices: 



ml :=cos 0 n !-sin ( 



2 



ml a 


4 

"a 


: ml a 2 ' n a 2 


4 ml a 2 n a 2 


"a 4 


ml a 4 


2 ml a 2 ' "a 2 


4 m 'a 2 "a 2 


m, a * "a 2 


ml a 2 ‘ n a 2 


"a' 


-4 ml a 2 - n a 2 


ml a 2 ' "a 2 


ml a 2 "a 2 


- 2 m * a 2 ' "a 2 


[ m *a 2— n a 2 ]~ 


ml a 3 "a 


-"''a' "a 3 


ml a "a 3 " ml a "'"a 


2 (ml a a a 3 - ml a 3 a a ] 


ml a "a 3 


- m *a 3 ' n a 


ml a 3n a- ml a' n a 3 


2 [ 3 n a- ml a n a 3 ] 



Ply Stiffness Calculation: 

Calculates the zero direction stiffness matrix for the various ply states (undamaged, matrix damaged, and 
fiber damaged) 



E xx q 


v yy q E xx q n 


'-'Vn, 


1 ~ v x.\q v y> q 


v a.\,' E » q 


E yv 

0 


1 — v v w 

NX q -q 


?, 

V 

> 

1 


0 


0 E cc 

ss q J 



Creates the stiffness matrices for every combination of ply angle and ply state. Qnew is a 3 x 4 matrix 



q i,i 



q i .2 



Q r °tq , a , Qrot q,a 3 Q ™q,a 5 
^ ro[ q,a 3 ^ rot q,a 2 ^ rot q,a 6 
Q rot q,a 5 Q rot q,a 6 Qrot q,a 4 



Undamaged Elastic Stiffness of Laminate: 

Q Laminate ;= N 0+ N 90 -f N n45 -fN ^5 ^O^new,, *h N 90' Q new, 2 +N n45 ' Q nevv, 3 -h N p45 Q new, 4 ] 



Q laminate 3 



17.575 5.256 -0.0381 

5 256 17.726 -0 038 

-0.038 -0.038 6.198 



°GPa S :-Q | am i nat e" 



P - 1 

c Eglass 1 - c 

1,1 
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Effective Thickness Calculations: 



A stiffener'^ stiffener E carbon 
1 eff ,_l Ld— — -> - 1 R . „ -p— ; — 
n — ^ ) K S h e n t Egi^s 



H efT := V R shell 1 eff 

- = L shell 



t e ff = 0.035Tn 
H ef f=0 I94*m 



t L = 0 0128*m 
H = 0 1 17*m 



, = 2 159 



N u m ber waves • -2 H e ff ' = — — 

s ett Number waves 2 H efr =0.21«m 

Flexural Stiffness: 

1. Indexing variables fr:=i.. Total _piys-t- 1 

ag :=1 Total_plys 

2. Establishing an array of thickness positions for each ply 

3. Assigns the appropriate stiffness matrix based on the lay-up sequence 



Q fiex 



gg 



Qnew, j if sequence g£ =0-deg 
Q new , 2 if sequence gg =90 deg 
Qnew| 3 if sequence ge =45 deg 
Q new i 4 otherwise 



Stress Space Strength Parameters: 

A separate F must be calculated for each ply state 



xx q xx pnmeq 






/XX -XX nnme YY YY 
y q primeq q 



primeq 



YY q YY primeq 



0 



0 




1 1 



XX q XX primeq 
1 _ 1 
YY q YV pnmeq 



/ 
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Strain Space Strength Parameters 

1 . Strain strength parameter for each ply state 





F ^ 1,1 Q q 1,1 Q q 1,2 


+ ~ F ^ 1,2 Q q 1,1 Q q 1,2 


+ F a ( Q • Q 

*q 1,2 L q u q 2,: 


+ F ^ 2 , 2 [ Q q ., 2] 2 


F ^ 2,2 Q q 1,2' Q q 2,2 


Fa q i,i ' Q q i.i' Q q 1,2 ' 


Fa q u '[ Q q 1.2] 2 


+ Fa q 1 ,2 [ Q q 1.1 Q q 2,2^"[ Q q 1,2] ] 


+2 ' Fa q 1,2' Q q 1.2 


^ F a q2 , : - V 2 Q q 2.2 


+ F -q 2,2 [ Q q 2,2] 


0 


0 



: w[ W 



Fb q , Q q ,.1 + Q q 1.2 
F b q i ^q 1,2^ F b q ; ^q 2,2 



2 Calculating G matrices for all combinations of ply state and angle: 





G ltrans :: 



Rotation Ij 


G baselj | 2 
G baselj ^ 3 J 

Gbase 'l l.l] 
G basel, , , 


Rotation l-> 






G basel | ( 

G basel 1 , , 

' Gtoel; , ,1 
G base 1 2 t 


Rotation U 






G basely t 
G basel 1 , , 
G basel 3 
G basel 3 


Rotation I-, 






G basel 3 j , 

G basel 3 , 

j, 3J 

G basel, , 

G basel | ^ 


Rotation ly 




G basel | | , 

G basel, , , 
j,jJ 

G basel 2 j j 
G basel 2 ^ 


Rotation 1^- 




G basel 2 , , 

G basel-) , , 

“ ->,3J 

G basel 3 ( J 
G basel 3 


Rotation 1 3 - 


G base 1 3 12 
G basel 3 , , 
G basel , L1 ‘ 
G basel , 


Rotation^ 


G basel , t , 
Q basel, 3 3 
G basel 2 ( j 
G basel 2 


Rotation^ 


G baseH j , 
G basel2 , 3 
G base 1 3 u ' 
G basel 3 t 


Rotation!^ 





Rotation^ 



*3 


G base2 , 


Rotation2j 


G base22 j 
G base2 : 2 


Rotation^ 


' G base2 3 
G base23 ? 


Rotation2 4 


’ G base2, ' 
G base2, 2 


Rotation2 4 


G base22 j 
G base22 2 


Rotation2 4 


G base23 j 
G base2 3 
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